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were  approximately  the  velocity  of  sound  in  water.  The  hyperbolic  signals  also  showed  dispersion  properties 
similar  to  acoustic  normal  modes  in  shallow  water.  These  observations  indicate  that  the  signals  responsible  for  the 
hyperbolic  reflections  propagate  as  normal  modes  within  the  water  layer,  with  very  limited  penetration  of  the  sea¬ 
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HYPERBOLIC  REFLECTIONS  ON 
BEAUFORT  SEA  SEISMIC  RECORDS 


K.G.  Neave,  P.V.  Sellmann  and  A.  Delaney 


INTRODUCTION 

Hyperbolic  reflections  are  a  common  feature  on 
marine  seismic  records  obtained  during  oil  exploration 
in  the  Beaufort  Sea  (Fig.  la).  These  reflections  were 
also  found  in  abundance  by  Reimnitz  et  al.  (1972) 
on  their  sub-bottom  seismic  profiles.  The  hyperbolic 
reflections  are  important  for  two  reasons:  1 )  they 
are  a  noise  problem,  and  2)  they  are  an  indication  of 
a  local  contrast  in  sediment  properties. 

This  noise  problem  on  oil  industry  records  that  are 
saturated  with  hyperbolas  would  require  expensive 
extra  record  processing  to  delineate  the  deeper  geo¬ 
logical  structures.  The  hyperbolas  that  saturate  the 
USGS  sub-bottom  profiles  taken  near  Prudhoe  Bay 
(Fig.  1  b)  also  contribute  to  the  noise  that  obscures 
any  structural  features  that  may  be  present  in  these 
shallower  sediments. 

These  reflections  are  also  important  since  they 
originate  from  features  within  the  sediments  that  are 
density  or  elastic  property  anomalies.  Reimnitz  et  al. 
(1972)  proposed  two  possible  sources  for  the  reflec¬ 
tions:  massive  ice  and  erratic  ooulders.  The  available 
evidence  suggested  to  them  that  massive  ice  was  the 
most  probable  cause.  Rogers  and  Morack  (1978) 
added  gas  pockets  and  Grantz  added  ice  gouges 
(Toimil  1978)  to  the  list  of  possible  sources.  Identifying 
the  source  of  these  anomalies  is  important  in  assessing 
hazards  to  construction  or  drilling  activities  in  this 
coastal  region. 

Velocity  measurements  for  the  materials  which 
surround  the  structural  features  in  question  can  be 
made  from  the  hyperbolic  reflections  since  their  shape 
is  a  function  of  velocity.  We  selected  a  curve-matching 
technique  to  analyze  the  type  of  records  available  for 


this  study.  A  projector  system  was  designed  which 
permitted  rapid  curve  matching  for  the  range  of 
hyperbolas  on  the  records. 

This  study  examines  the  velocities  and  dispersion 
characteristics  of  the  hyperbolic  reflection  signals 
and  proposes  that  many  originate  from  anomalies 
within  several  meters  of  the  seabed.  Their  common 
occurrence  suggests  that  they  must  originate  from 
abundant  features  whose  material  properties  contrast 
significantly  with  those  of  the  surrounding  material. 
The  anomalous  reflections  originate  from  both  point 
and  linear  sources.  The  low  velocities  suggest  that 
these  reflections  are  not  likely  to  be  structures  incor¬ 
porated  in  ice-bonded  sediment. 

We  propose  that  these  reflections  originate  from 
the  filled  or  unfilled  common  gouges  and  troughs 
formed  in  the  seabed  on  much  of  the  continental 
shelf  by  impingement  of  floating  ice.  These  gouges 
commonly  occur  in  dense,  overconsolidated,  fine¬ 
grained  sediment  and  can  be  filled  with  lower  density 
material  (Reimnitz  and  Barnes  1974,  Reimnitz  and 
Toimil  1 979).  The  filling  can  occur  rapidly  in  the 
course  of  periodic  storms  during  unusual,  ice-free 
conditions.  Seabed  sediments  are  redistributed  at 
this  time  by  increased  wave  and  current  action  on  the 
shelf.  Barnes  and  Reimnitz  (1979)  observed  that  one 
such  storm  caused  extensive  obliteration  of  ice  gouges 
in  water  at  least  1 3  m  deep  in  their  study  area  and 
caused  sediment  ponding  and  gouge  filling  at  even 
greater  depths.  The  shape  and  size  of  the  gouges,  the 
potential  existence  of  buried  gouge  structures,  and  the 
variations  in  bed  sediment  properties  caused  by  re¬ 
working  of  the  seabed  by  the  ice  all  fall  within  the 
range  of  seabed  modifications  which  can  reflect  energy 
and  cause  hyperbolas. 


METHODS  OF  ANALYSIS 


Marine  seismic  records 

Marine  seismic  records  were  analyzed  to  determine 
the  properties  of  the  hyperbolic  reflectors.  The  oil 
industry  data  we  used  were  recorded  by  Geophysical 
Service  Inc.  as  part  of  their  offshore  oil  and  gas  ex¬ 
plorations.  The  energy  source  was  a  tuned  air  gun 
array  with  approximately  I  5,000  cm5  (900  in.* ) 
capacity.  The  air  guns  were  towed  3  m  below  the 
surface  while  the  receiver  cable  was  at  a  depth  of  4.5  m. 
Thirty  hydrophones  spaced  1 .2  m  apart  made  up  a 
receiver  group  which  fed  into  a  single  channel  on  the 
recorder.  Forty-eight  groups  spaced  50  m  apart  made 
up  the  complete  2.35-km-long  cable.  For  recording 
the  data,  they  used  a  DFS  IV  digital  system  with  a 
4-ms  sampling  rate.  The  filter  settings  were  8  to  62  Hz. 

In  this  study,  only  the  early  arrivals  were  analyzed; 
therefore,  data  processing  on  playback  was  rudimentary. 
Two-second-long  monitor  records  were  printed  in  the 
variable  area  format  with  increased  early  gain  settings. 
There  was  no  mixing,  stacking,  or  gathering  of  the 
traces  during  playback.  There  was  adequate  coverage 
for  a  reconnaissance  survey,  with  a  record  for  each 
kilometer.  Since  the  shot  point  interval  was  50  m, 
we  only  required  the  playback  of  every  20th  shot 
point. 

The  signature  of  the  source  pulse  from  the  air  guns 
and  the  signal  modification  by  the  recording  system 
can  be  seen  in  the  recordings  of  signals  which  have  a 
simple  travel  path  through  the  water.  A  reproduction 
of  the  direct  wave  through  water  (Fig.  2a)  shows  that 
a  complex  waveform  is  produced  by  the  system; 
however,  the  complexity  can  be  attributed  mainly  to 
the  geometric  filtering  of  the  source  and  receiver.  In 
Figure  2b  the  reproduction  of  the  sea  floor  reflection 
in  deep  water  has  a  short  and  simple  signature.  The 
waveform  in  Figure  2b  should  be  the  appropriate 
form  for  signals  that  are  transmitted  and  received 
along  rays  that  are  nearly  perpendicular  to  the  array. 

The  oil  exploration  records  contain  strong  signals 
which  were  identified  as  normal  modes  of  trapped 
acoustic  energy  in  the  water  layer.  We  propose  that 
the  anomalous  reflections  of  normal  mode  energy  are 
involved  in  the  generation  of  the  hyperbolic  signals, 
and  we  have  made  comparisons  of  the  properties  of 
the  normal  modes-velocity,  frequency  content  and 
skin  depth— with  the  hyperbolic  reflections  to  test 
this  hypothesis. 

C.L.  Pekeris  established  a  theory  for  the  propagation 
of  the  acoustic  normal  modes  in  water  (Ewing  et  al. 

1 957,  p.  1 26).  His  calculations  show  that  the  largest 
amplitude  interval  of  these  signals  is  the  Airy  phase 
of  the  fundamental  mode.  His  dispersion  curves  were 
useful  for  predicting  the  velocity  and  frequency  to  be 
expected  on  the  records.  For  a  sea  floor  sediment 
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velocity  of  1 .8  km/s,  the  Airy  phase  of  the  funda¬ 
mental  mode  travels  at  a  group  velocity  that  is  6% 
less  than  the  water  velocity,  with  a  frequency  fj 
(in  Hz)  that  is  inversely  proportional  to  the  water 
depth  H  (in  meters): 

f ,  =  0.77/7/.  (1) 

Since  the  water  depth  is  known,  it  is  possible  to  pre¬ 
dict  this  frequency. 

Frequency  measurements  were  made  to  verify  the 
theoretical  analysis  for  the  normal  modes  and  to 
confirm  the  Airy  phase  identification  for  both  the 
water  wave  and  the  hyperbolic  reflections.  First,  the 
signal  period  was  determined  by  measuring  the  time 
interval  between  consecutive  peaks  on  the  record. 
Then,  the  dominant  frequency  was  calculated  from 
the  signal  period.  It  can  be  shown  that  this  measure¬ 
ment  method  gives  the  frequency  at  which  the 
Fourier  transform  of  the  signal  has  the  maximum 
amplitude.  For  a  long  burst  of  sinusoidal  energy, 
such  as  is  encountered  for  the  Airy  phase,  there 
would  be  a  narrow  peak  in  the  frequency  spectrum. 
Accordingly,  the  Airy  phase  signal  can  be  described 
approximately  by  the  frequency  at  which  the  fre¬ 
quency  spectrum  peaks. 

The  normal  modes  propagate  as  multiple  reflec¬ 
tions  within  the  water  layer  along  with  matching 
disturbances  within  the  seabed  sediments.  Ewing 
et  al.  (1957,  p.  79)  calculated  that  the  amplitude  A 
of  these  disturbances  suffered  an  inverse  exponential 
decay  with  depth  Z  below  the  sea  floor,  such  that 

A=40exp(-Z/Z0).  (2) 

The  quantity  ZQ  is  a  characteristic  attenuation  depth 
called  the  skin  depth.  It  is  related  to  the  frequency  of 
the  signal  f  by  the  expression 

Z0  =  \clf\\}-(c!V2)2)-'/2  (3) 

where  c  is  the  horizontal  phase  velocity  and  V  2  is 
the  compressional  wave  velocity  in  the  sediments. 

The  normal  modes  may  be  viewed  as  being  carried  in 
a  waveguide  that  is  made  up  of  the  water  layer  plus 
the  sediments  near  the  sea  floor  down  to  the  skin 
depth.  Undisturbed  transmission  of  a  normal  mode 
would  require  a  homogeneous  layer  of  sediments 
with  a  thickness  of  the  skin  depth  ZQ.  Disturbances 
or  unexpected  reflections  of  the  normal  modes  would 
be  created  by  any  significant  irregularity  in  the  seabed 
or  an  anomaly  in  the  clastic  properties  of  the  sedi¬ 
ments  within  the  skin  depth.  Skin  depths  for  the 
Airy  phases  were  determined  using  eq  3  for  compar¬ 
ison  with  gouge  depths. 


a.  The  direct  wave  in  water  recorded  by  the  four  hydrophone  groups  closest 
to  the  source.  This  shows  the  system  characteristics  /'  '  signals  traveling 
parallel  to  the  array. 
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1  ..  -  t»  —  p  -120m  b.  Sea  floor  reflection  in  deep  water  recorded  by  the  same  four  hydrophones. 

I  |  |  This  shows  the  system  response  for  a  wave  traveling  nearly  perpendicular  to 
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Figure  2.  Source  and  recording  system  response  for  the  oil  industry  records. 


Velocity  measurements  are  a  basic  tool  for  deter¬ 
mining  the  identity  of  the  signals  and  analyzing  the 
propagation  path.  For  the  water  wave  on  the  oil  in¬ 
dustry  records,  the  velocity  of  the  peak  amplitude  of 
the  modes  was  determined  using  a  standard  time- 
distance  analysis.  Flowever,  velocity  determinations 
for  the  hyperbolas  were  not  so  straightforward. 
Velocities  can  be  determined  from  the  hyperbolic 
reflections  by  a  curve-matching  technique,  using  the 
projector  described  in  the  Appendix.  A  known  curve 
can  be  found  for  each  hyperbola  on  the  records  and 
the  parameters  of  the  known  curve  will  convert  to  a 
velocity. 

In  order  to  make  the  measurements  from  the  hy¬ 
perbolas,  an  assumption  was  required  about  the  re¬ 
flector  geometry,  since  the  velocity  is  partly  controlled 
by  what  shape  is  assumed.  Flyperbolas  could  originate 
from  any  reflector  shape  between  a  point  anomaly  and 
a  plane  reflector.  The  reflection  sources  were  assumed 
to  be  linear  features  in  the  seabed  sediments  in  accor¬ 
dance  with  the  proposed  reflection  mechanism. 

The  oil  exploration  recording  system  has  a  signal 
source  followed  by  a  linear  array  of  receivers  (Fig. 

3a).  The  reflecting  anomaly  is  depicted  as  a  line 
which  can  have  an  arbitrary  orientation  on  the  bot¬ 
tom  with  respect  to  the  receiver  array.  For  shallow 
water,  the  vertical  component  of  the  propagation 
path  can  be  ignored;  thus  the  problem  is  two-dimen¬ 
sional.  In  this  diagram,  the  real  signal  source  is 
labeled  S  and  the  image  source  behind  the  reflector 
is  labeled  S'.  The  X-axis  lies  along  the  line  of  receiv¬ 
ers,  and  the  origin  O  is  placed  at  the  point  on  the 
axis  closest  to  the  image  source  S’.  The  equation 
that  describes  the  propagation  time  t  to  a  receiver  at 
a  position  s  is 

t  =  (1/l/0)  \/x2  +  (2 y)2  (4) 


where  VQ  is  the  propagation  velocity  and  2 y  is  the 
distance  OS' .  On  a  seismogram,  the  plotting  axes  are 
distance  £  and  time  t.  The  coordinate  £  on  the  record 
is  related  to  the  real  distance  coordinate  of  Figure  10a 
by  a  known  scaling  factor  k : 

£  =  6*.  (5) 

Combining  eq  4  and  5  gives  the  equation  for  reflec¬ 
tion  times  on  a  seismogram: 


This  is  the  same  type  of  equation  as  in  the  Appendix 
for  the  set  of  hyperbolas  which  are  available  for  curve 
matching: 

(t/a.y-di/b;)2-  1  (7) 

where  a-t  and  bt  are  constants  which  are  determined 
by  the  projector  system  and  /  =  1  to  n.  If  projector 
curve /  matches  the  curve  on  the  seismogram,  the 
parameters  from  eq  6  and  7  must  be  equal  for  /  =/. 
Therefore,  there  is  a  straightforward  solution  for  the 
distance y  and  the  velocity  V0: 


y=(  M2k)b., 

(8) 

V'0  =  OM)(*j/«j). 

(9) 

Equation  9  was  used  to  obtain  velocities  from  the  oil 
exploration  records. 

We  expect  one  of  the  assumptions  used  in  the 
velocity  analysis  to  introduce  observational  errors  in 
the  hyperbola  velocities.  A  random  scatter  in  the 
velocity  readings  will  result  if  the  reflecting  anomalies 
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a.  4  ''Dlui  seabed  reflector  A-A '  near  the  oil  explora¬ 
tion  recording  system.  O  is  the  origin,  S  is  the  source, 
S'  is  the  image  source  and  R  is  one  of  the  receivers. 


b.  A  point  reflector  in  the  seabed  near  the  sub-bottom 
profiling  system.  The  symbols  are  the  same  as  in  a, 
with  P  being  the  point  reflector. 


Figure  3.  Plan  view  of  reflection  geometry  for  marine  seismic  records. 


are  curved  rather  than  straight.  Reflections  from 
anomalies  that  are  concave  toward  the  receive,  airay 
will  give  an  exaggerated  velocity,  while  reflections 
from  convex  anomalies  will  yield  an  underestimate 
of  the  velocity.  A  random  sample  of  anomaly  shapes 
will  therefore  yield  a  corresponding  scatter  in  velocity 
estimates;  however,  the  average  of  a  number  of  mea¬ 
surements  will  lead  to  a  good  estimate  of  the  propa¬ 
gation  velocity. 

Several  other  aspects  of  the  hyperbolic  reflectors 
on  oil  industry  records  can  be  seen  on  a  processed 
section.  The  hyperbolas  are  strong  signals,  frequently 
with  larger  amplitudes  than  the  sedimentary  horizon 
reflections.  They  are  similar  to  the  sedimentary  re¬ 
flections  and  we  suggest  that,  because  of  their  large 
amplitude,  they  will  not  be  removed  by  common- 
depth-point  gathers  with  normal- move-out  corrections. 
In  this  survey,  there  is  24-fold,  common-depth-point 
coverage.  Each  trace  on  the  section  is  made  by  add¬ 
ing  single  traces  from  24  adjacent  shot  points  with 
normal-move-out  corrections  to  compensate  for  the 
hyperbolic  shape  of  the  reflections.  The  comoensat- 
ing  times  are  calculated  using  an  appropriate  velocity 
for  the  deeper  sediments:  2.0  to  4.5  km/s;  however, 
it  is  possible  for  a  large  amplitude  hyperbola  with  a 
different  velocity  to  appear  on  the  traces  after  stack¬ 
ing.  When  the  traces  are  printed  together  as  a  section, 
the  hyperbolic  reflectors  which  have  not  been  sup¬ 
pressed  completely  by  the  processing  will  appear  as 
dipping  and  intersecting  features  along  with  the  sedi¬ 
mentary  structures.  The  sections  provide  continuous 
coverage  along  the  survey  line  so  that  the  length  and 
orientation  of  the  hyperbolic  reflectors  can  be  fol¬ 
lowed  over  considerable  distances.  This  is  a  great 


improvement  over  the  individual  monitor  records, 
which  can  cover  only  1.2  km  of  any  reflector.  A 
representative  sample  of  seismic  sections  from  the 
continental  shelf  west  of  Prudhoe  Bay  was  examined 
to  identify  the  hyperbolic  reflectors  and  find  the 
length  of  the  anomalies. 

Seismic  sub-bottom  profiles 

Hyperbola  data  were  also  collected  from  sub-bottom 
profiles  made  by  Reimnitz  et  al.  (1972).  Their  system 
consisted  of  an  arc  source  (sparker)  and  a  Gifft  fac¬ 
simile  recorder.  The  sparker  releases  500-)  impulses 
into  the  water  at  K-s  intervals.  The  receiver  used  elec¬ 
tronic  filters  with  a  passband  between  430  and  960 
Hz  and  printed  analog  records  'A  s  long. 

An  image  of  the  source  pulse  can  be  seen  in  the 
bottom  reflection  on  a  deep  water  record  (Fig.  4). 

The  reflection  energy  has  a  duration  of  6  ms  and  a 
dominant  frequency  of  approximately  400  Hz  (deter¬ 
mined  by  measuring  the  time  between  peaks  on  the 
record). 

The  sub-bottom  profiles  are  shot  with  a  single  re¬ 
ceiver  instead  of  an  array  as  used  for  the  oil  exploration 
records.  In  contrast  to  the  oil  industry  data,  there  arc 
no  problems  in  determining  the  shape  of  the  reflectors 
on  the  sub-bottom  profiles.  The  source  and  receiver 
have  a  common  location  and,  as  a  result,  the  outgoing 
and  returning  rays  have  the  same  path.  For  this  rea¬ 
son,  a  point  reflector  creates  a  hyperbola  on  the  record, 
whereas  a  finite  plane  reflector  shows  up  as  a  linear 
feature  with  half-hyperbolas  on  each  end. 

Figure  3b  shows  the  sub-bottom  profiler  system, 
with  the  reflecting  anomaly  depicted  as  a  point  re¬ 
flector.  The  source  and  receiver  are  at  the  point  S. 
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Figure  4.  Part  of  a  sub-bottom  profile  taken  north  of  Oiiktok  Point. 
The  water  depth  is  20  m. 


The  X-axis  lies  along  the  ship's  track  and  the  origin 
O  is  placed  where  the  X-axis  is  closest  to  the  point 
reflector  P.  To  describe  the  travel  time  t  for  the  re¬ 
flected  signal,  the  equation  is 

t~-{7iv0)sj^r-2  (io) 

where  x  is  the  distance  05,  and  i  is  the  distance  OP. 
When  a  profile  is  recorded,  the  two  axes  on  the  plot 
arc  distance  |  and  time  t.  The  coordinate  £  on  the 
record  is  related  to  the  real  distance  coordinate  x  by 
two  constants,  the  ship’s  speed  r  and  the  paper  speed 
Q- 

£=(«?/'')*•  (11) 

Combining  eq  10  and  1 1  gives  the  equation  for  re¬ 
flection  times  on  the  sub-bottom  profiles: 

<'2> 

This  equation  is  similar  to  the  set  of  curve-matching 
equations  in  the  Appendix.  The  same  technique  as 
before  is  used  to  solve  for  the  two  unknowns,  a  and 
^o- 

z=(r/q)bm  (13) 

V0--2{r/g)(bm/aJ.  (14) 

equation  14  was  used  to  find  the  velocities  for  the 
reflections  on  the  sub-bottom  profiles. 

A  comparison  between  the  dispersion  properties 


of  the  hyperbolic  reflections  and  those  of  the  normal 
modes  in  the  water  layer  was  made  using  a  synthetic 
seismogram  generated  by  the  hyperbola  projector  sys¬ 
tem.  The  signal  frequency  and  velocities  were  taken 
from  Pekeris’  theory  (Ewing  etal.  1957).  A  water 
depth  of  3  m  was  chosen,  and  the  maximum  ampli¬ 
tude  part  of  the  signal,  the  Airy  phase  of  the  funda¬ 
mental  mode,  was  used  for  the  plot.  The  resulting 
record,  showing  the  noticeable  influence  of  dispersion, 
can  be  seen  in  Figure  5.  Dispersion  caused  the  peaks 
and  troughs  of  the  reflections  to  cut  across  the  wave 
packet  at  an  angle,  producing  a  pattern  known  as 
shingling.  Even  though  the  signal  energy  runs  contin¬ 
uously  down  the  record,  each  individual  peak  or 
trough  can  only  be  traced  for  about  40  ms  before  it 
is  replaced  by  a  later  one.  The  hyperbolas  were  exam¬ 
ined  for  evidence  of  dispersion  because  it  is  a  charac¬ 
teristic  feature  of  normal  modes  which  is  not  observed 
for  ordinary  congressional  wave  reflections  in  the 
sub-bottom  sediments. 

Multiple  reflections  can  be  used  as  an  indication  of 
the  reflection  geometry  on  the  profile  records.  Mul¬ 
tiples  are  generated  as  the  signal  bounces  between  any 
strong  reflector  at  depth  and  the  strongly  reflecting 
air/water  interface.  Typical  multiples  are  visible  on 
the  record  in  Figure  4.  They  arrive  at  twice  and  three 
times  the  reflection  time  of  the  primary  reflection 
from  the  seabed.  If  there  is  a  strong  point  reflector 
at  depth  which  generates  a  hyperbolic  reflection,  it  is 
then  possible  to  get  a  multiple  reflection;  however,  the 
arrival  time  for  the  multiple  will  be  less  than  twice 
the  primary  time  when  the  point  reflector  is  not  di¬ 
rectly  under  the  ship’s  track.  Therefore,  we  searched 
the  records  to  find  hyperbolic  multiples  at  times 


Figure  5.  Synthetic  seismogram  for  a  normal  mode  in  water  3  m  deep  reflected  by  a 
point  reflector.  Dispersion  causes  the  peaks  and  troughs  to  cut  diagonally  across  the 
wave  packet. 


intermediate  between  the  primary  arrival  and  twice 
the  primary  time. 

RESULTS  AND  DISCUSSION 

Distribution  of  hyperbolic  reflections 

Hyperbolic  reflections  are  found  on  nearly  all  the 
oil  exploration  lines  shown  in  Figure  6.  They  occur 
across  the  entire  survey  area  from  Prudhoe  Bay  to 
Demarcation  Bay.  Additional  data  taken  from  the 


Prudhoe  Bay  to  Pt.  Barrow  area  also  show  places  with 
abundant  hyperbolas.  These  areas  are  sporadically 
distributed  between  the  continental  shelf  break  (near 
the  100-m  isobath)  and  the  shallowest  parts  of  the 
seismic  lines  (usually  in  about  15  m  of  water),  with 
slightly  increased  frequency  near  the  coastline. 

Reimnitz  et  al.  (1972)  found  hyperbolic  reflections 
on  the  sub-bottom  profiles  over  a  more  restricted  range, 
as  shown  on  their  map  in  Figure  7.  Their  distribution 
was  confined  to  the  shallower  water  near  the  shore  of 
the  mainland  and  around  the  barrier  islands.  They 


Figure  6.  Distribution  of  hyperbolic  reflections  on  oil  industry  seismic  data.  The  seismic  lines  are  plotted  as  narrow 
lines  and  the  locations  of  hyperbolic  reflections  are  identified  by  wide  lines. 
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Distribution  of  hyperbolic  reflections  on  the  sub-bottom  profiles  of  Reimnitz  et  at.  (1972). 


observed  few  hyperbolic  reflections  where  the  water 
is  deeper  than  6  m.  This  apparent  dependence  on 
water  depth  could  be  related  to  the  frequency  filters 
of  their  receiver  system.  Their  filtering  discriminated 
against  low  frequencies  which  are  carried  by  the  nor¬ 
mal  modes  in  deeper  water.  The  effect  of  filtering  is 
discussed  in  a  following  section. 

Hyperbolas  on  oil  exploration  records 

A  direct  comparison  can  be  made  between  the 
water  wave  normal  modes  and  the  hyperbola  signals 
in  Figure  la  (the  oil  exploration  records).  There  is 
noticeable  resemblance  between  their  frequency  coti: 
tent  and  pulse  length.  Figure  8  shows  reproductions 
of  both  types  of  signal  so  that  their  similarity  can  be 
observed.  They  both  have  approximately  0.1  -s  sinu¬ 
soidal  bursts  of  energy  with  similar  frequency  content. 

The  dominant  frequency  calculated  for  the  water 
wave  modes  from  measurements  made  on  some  of 
the  data  from  the  western  part  of  the  survey  area  was 
about  50  Hz,  compared  with  45  Hz  for  the  hyperbolic 
data  (Fig.  9).  This  is  only  a  10%  difference,  which 
might  be  accounted  for  by  interference  and  filtering 
during  the  reflection  process. 

The  depth  of  signal  penetration  into  the  seabed 
can  be  calculated  using  the  frequency  value  of  45  Hz. 
Based  on  eq  3,  the  normal  modes  will  have  a  skin 
depth  of  1 2  m.  This  is  the  same  order  of  magnitude 
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f  iqure  8.  Water  wave  normal  modes  compared  to  a  hyper¬ 
bolic  reflection.  Traces  A  and  B  (from  the  sample  record 
of  Figure  la)  have  maximum  amplitudes  for  the  Airy  phase 
of  the  fundamental  mode.  The  group  of  traces  labeled  C 
contain  part  of  a  hyperbolic  reflection  from  the  same  re¬ 
cord  which  appears  similar  to  the  A  iry  phase  in  traces  A 
and  B. 


as  the  penetration  depth  of  the  deeper  gouges  reported 
for  the  Beaufort  Sea:  about  5.5  to  6.5  m  (Barnes  and 
Hopkins  1978). 

The  calculated  hyperbolic  velocities  (Fig.  10)  gave 
an  average  of  1 .3  km/s.  This  average  velocity  is  within 
the  experimental  error  of  the  unreflected  normal  mode 
velocity:  1.35  km/s.  The  range  in  the  velocities  of  0.8 
to  1 .9  km/s  can  be  attributed  to  curvature  of  the  re¬ 
flecting  anomalies. 

It  appears  that  the  noise  problem  created  by  hyper¬ 
bolas  on  the  marine  records  can  be  seen  in  the  pro¬ 
cessed  seismic  section  reproduced  in  Figure  1 1.  Only 
pne  sedimentary  horizon  is  continuously  visible  across 
the  section  at  1 .6  s.  Several  parallel  reflectors  can  be 
found  intermittently  between  0.7  and  1.6  s.  All  the 
other  intersecting,  dipping  and  curving  lineaments  on 
the  section  appear  to  result  from  hyperbolic  reflections 
which  were  not  adequately  suppressed  by  the  process¬ 
ing  techniques. 

The  seismic  section  in  Figure  1 1  can  be  treated  as 
a  map  of  ice  gouges  if  the  sedimentary  structures  are 
ignored.  A  high  gouge  density  is  apparent  for  this 
sample  record  which  was  shot  in  35  m  of  water  off 
Cape  Simpson.  Continuous  gouges  up  to  4  km  long 
are  present.  Their  dominant  frequency  of  20  Hz  im¬ 
plies  that  the  skin  depth  of  the  signal  is  27  m. 


Wot»r  Wav* 
(normal  mode) 
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Figure  9.  Dominant  frequency  measure¬ 
ments  of  signals  on  the  oil  exploration 
records. 
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Figure  I  I.  A  processed  marine  seismic  section  shot  by  GS/  near  Cape  Simpson,  Alaska,  and  obtained  from  the  U.S. 
Government.  Some  horizontal  sedimentary  structure  is  visible  between  0. 7  and  1.8  s.  The  dipping  and  intersecting 
lineations  are  interpreted  as  ice  gouges  in  the  seabed.  (A  nonpropi  ietary  section  released  js  part  ot  the  seismic  data 
generated  in  the  NPRA.) 
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Figure  1 2.  Hyperbola  measurements  from  sub-bottom  profiles  of  Prudhoe  Bay  and 
Simpson  Lagoon. 
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Hyperbolas  on  sub-bottom  profiles 

Most  of  the  hyperbolic  reflections  on  the  sub¬ 
bottom  profiles  appear  to  be  point  reflections.  In  a 
few  cases,  linear  features  were  found  which  indicate 
anomalies  up  to  60  m  long;  however,  only  the  point 
reflectors  were  used  for  the  velocity  analysis.  Two 
lines  from  Prudhoe  Bay  and  one  line  from  Simpson 
Lagoon  were  used  for  the  velocity  measurements. 

An  average  velocity  of  1 .5  km/s  was  observed,  with 
values  ranging  from  l  .l  to  1.9  km/s  (Fig.  12a).  The 
scatter  in  the  readings  can  be  attributed  to  the  reso¬ 
lution  of  the  measuring  system.  A  possible  error  of 
i_0.4  km/s  was  determined  for  these  measurements 
in  the  Appendix.  The  average  velocity  from  the  hyper 
bolas  is  close  to  the  water  wave  velocity,  1 .4  km/s,  and 
also  close  to  the  average  velocity  in  the  sub-bottom 
sediments  determined  by  Rogers  and  Morack  (1978): 

1 .8  km/s.  The  resolution  of  the  observations  from  the 
hyperbolic  reflections  is  not  great  enough  to  exclude 
consideration  of  the  sediment  velocity. 

The  distance  from  the  source  of  the  reflectors  fell 
between  40  and  160  m  (Fig.  1 2b).  The  lower  limit 
was  determined  by  the  projector  measurement  method 
and  not  by  the  reflector  distribution. 

Shingling  of  the  hyperbolas  is  clearly  evident  on  the 
sub-bottom  profiles.  Figure  1 3  shows  examples  where 
each  peak  (represented  by  one  line)  can  be  traced  for 
20  to  40  ms,  whereas  the  energy  extends  over  approx¬ 
imately  75  ms.  The  peaks  propagate  at  a  higher  veloci¬ 
ty  than  the  energy  envelope,  as  expected  for  dispersion. 
This  is  the  same  type  of  behavior  predicted  in  the  syn¬ 
thetic  seismogram  for  a  reflected  normal  mode  (Fig. 

5).  The  dispersion  of  the  hyperbolas  is  contrary  to  the 
expected  signal  from  simple  concessional  waves  re¬ 
flecting  off  small  anomalies  buried  deep  in  the  sedi¬ 
ments. 

We  measured  the  dominant  frequency  of  the  hyper¬ 
bolas  to  determine  if  it  was  compatible  with  normal 
mode  propagation  in  the  water  layer.  The  average 


water  depth  for  the  chosen  survey  lines  was  2.5  m 
which  would  support  the  fundamental  mode  peak 
amplitudes  at  280  Hz.  The  histogram  of  the  dominant 
frequencies  of  the  hyperbolas  from  the  Prudhoe  Bay 
lines  (Fig  14)  shows  that  a  frequency  of  330  H/  is  most 
common.  This  illustrates  that  the  dominant  frequencies 
of  the  hyperbolic  reflections  are  consistent  with  normal 
modes  in  the  water  layer  and  somewhat  lower  than  the 
dominant  frequency  of  the  source:  400  Hr.  Agreement 
between  the  modal  theory  and  the  measurements  does 
not  need  to  be  exact  since  the  low-cut  filters  of  the 
receiver  were  set  at  430  Hr.  Because  of  this,  all  the 
low  frequency  components  of  the  signals  were  i  -duced 
in  amplitude,  and  the  hyperbolas  from  deeper  water 
were  suppressed  or  removed. 

The  common  frequency  of  330  Hr  for  the  Prudhoe 
Bay  reflections  can  be  used  to  estimate  the  depth  to 
which  the  normal  mode  signal  penetrated  the  sediments. 
In  this  case,  the  skin  depth  was  approximately  1 .5  m, 
a  depth  that  seems  to  be  compatible  with  the  smaller 
scale  of  gouging  in  shallow  water. 

No  significant  evidence  was  found  for  multiple 
reflections  from  the  hyperbolas.  For  50  of  the  most 
prominent  reflections,  there  were  only  ten  cases  where 
a  second  similar  signal  was  observed  between  the 
arrival  time  of  the  primary  and  twice  that  time.  This 
proportion  is  small  enough  that  we  interpret  it  as  a 
coincidence  of  primary  hyperbolas  rather  than  multiple 
reflections.  The  lack  of  multiple  reflections  from  these 
prominent  reflections  is  in  agreement  with  the  disper¬ 
sion  observation;  i.*.,  the  hyperbolic  signals  travel  in  a 
near-surface  waveguide  and  not  in  the  deeper  subsea 
sediments. 


CONCLUSIONS 

We  recognized  the  energy  carrier  for  the  hyperbolic 
reflections  on  the  oil  industry  records  as  acoustic 
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a.  A  sample  record. 
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b.  A  tracing  of  the  sample  record  to  emphasize  the  shingling, 
figure  13.  Hyperbolas  from  Simpson  Lagoon  showing  dispersion. 


Figure  14.  Dominant  frequency  of  hyper¬ 
bolas  on  sub-bottom  profiles  from  Prudhoe 
Bay  and  Simpson  Lagoon.  The  expected 
theoretical  normal-mode  periods  are  indi¬ 
cated,  as  well  as  the  filter  range  of  the  re¬ 
ceiver. 


normal  modes  in  the  water  layer.  We  suggest  that  the 
energy  is  most  likely  reflected  by  gouges  in  the  seabed 
caused  by  grounded  ice.  Gouges  filled  with  water  or 
soft,  recent  sediment  could  provide  enough  contrast 
in  material  properties  to  cause  a  reflection.  Identifi¬ 
cation  of  normal  mode  propagation  was  based  on  the 
similarity  of  the  signal  frequency  and  pulse  length  for 
the  observed  normal  modes  and  hyperbolic  reflections. 
Data  supporting  the  idea  of  the  gouges  being  the 
source  of  the  reflections  also  come  from  the  agree¬ 
ment  between  the  zones  where  gouges  are  anticipated 
and  the  observed  distribution  of  the  reflectors.  Fur¬ 
thermore,  the  skin  depth  of  12  m,  determined  from 
signal  frequencies,  would  accommodate  structures 
formed  by  the  deepest  gouges  observed,  as  well  as 
gouge  structures  that  may  have  been  incorporated  in 
the  bed  by  sedimentation.  Modification  of  the  sedi¬ 
ments  below  a  gouge  structure  could  also  act  to  in¬ 
crease  the  apparent  depth  of  the  structure.  Additional 
evidence  supporting  the  gouge  concept  is  that  an 
assumed  linear  geometry  for  the  reflectors  yields  a 
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reflection  velocity  which  agrees  with  normal  mode 
propagation  in  the  water  layer. 

Other  sources  for  the  reflections  (listed  in  the  intro¬ 
duction)  do  not  appear  to  provide  satisfactory  expla¬ 
nations  for  the  observations  made  as  part  of  this  study. 
The  observed  velocities  were  below  those  associated 
with  ice-bonded  marine  sediments,  so  ice  lenses  in 
frozen  sediment  do  not  provide  a  satisfactory  expla¬ 
nation.  The  observed  signal  frequencies  are  hard  to 
explain  if  gas  or  boulders  are  the  reflectors.  Com- 
pressional  wave  reflections  from  these  types  of  anom¬ 
alies  should  have  a  signal  with  a  frequency  content 
that  is  determined  by  the  size  of  the  anomaly  and 
that  should  not  necessarily  correspond  to  the  fre¬ 
quency  of  the  normal  modes.  In  addition,  the  velocity 
measurements  from  the  hyperbc'as  would  require 
calculations  based  on  a  point-source  geometry.  The 
linear  geometry  would  not  be  valid  for  the  oil  industry 
records.  As  a  result,  the  agreement  between  the  water- 
wave  velocity  and  the  hyperbola  velocities  would  have 
to  be  discounted. 

Analysis  of  the  sub-bottom  profiles  also  strongly 
points  toward  acoustic  normal  modes  in  the  water 
layer  as  the  propagation  mechanism  for  the  hyper¬ 
bolic  signals.  This  mechanism  provides  a  satisfactory 
explanation  for  the  dispersion  and  the  low  velocities 
that  were  measured.  When  the  signal  frequency  was 
compared  with  the  theoretical  normal  mode  frequency, 
the  agreement  was  good,  considering  the  type  of 
receiver  filtering  used.  Propagation  as  normal  modes, 
as  in  the  case  of  the  marine  records,  suggests  that  the 
anomalies  occur  very  near  the  seabed.  The  normal 
mode  signal  would  only  penetrate  'vl.Sm  into  the 
bed;  therefore,  it  does  not  seem  that  ice-bonded  perma¬ 
frost  anomalies  or  gas  concentrations  can  be  present 
at  this  depth  to  act  as  reflectors.  Boulders  lying  on 
or  just  below  the  seabed  could  act  as  point  reflectors 
for  the  normal  modes;  however,  the  boulders  could 
not  cause  the  linear  anomalies  which  appear  on  some 
records. 

There  are  two  practical  benefits  which  emerge  from 
the  gouge  explanation  of  the  hyperbolas.  First,  it 
may  be  possible  to  improve  the  signal-to-noise  ratio 
on  the  records.  This  could  be  accomplished  if  a  system 
could  be  designed  to  minimize  the  generation  and  re¬ 
ception  of  the  water  wave  normal  modes. 

The  second  potential  benefit  from  the  gouge  con¬ 
cept  is  the  possibility  of  locating  and  mapping  the 
major  gouge  tracks  on  the  continental  shelf  by  analysis 
of  existing  oil  exploration  records.  Similarly,  the 
smaller  gouges  and  associated  sedimentary  structures 
near  the  coast  could  be  mapped  using  the  sub-bottom 
profiler.  Even  relic  gouge  tracks  which  have  been 
filled  with  fresh  sediments  should  show  up  in  such  a 
study.  In  view  of  the  problems  that  ice  gouging  can 


present  to  development  along  the  Beaufort  Sea  coast, 
it  may  be  worthwhile  tQ  add  these  methods  to  the  list 
of  techniques  that  might  be  used  for  their  study. 
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APPENDIX  A:  HYPERBOLA  PROJECTOR 


The  projector  system  is  designed  to  generate  a  set 
of  standard  hyperbolas  tor  curve  matching  on  the 
records.  This  is  accomplished  by  shining  a  point 
source  ot  light  on  a  clear  sheet  which  has  been  etched 
with  a  set  of  ellipses  (Fig.  Al).  The  method  requires 
the  light  source  to  be  mounted  at  the  origin  0  of 
f  igure  A2  at  a  height  H  above  the  record,  and  the 
set  of  ellipses  suspended  with  their  centers  on  the 
r-axis  at  t  -  r.  The  set  of  ellipses  can  be  described 
in  this  coordinate  system  as 

U/d,)2  +  (Hla,)2  -  I  (Al) 

where  and  d,  are  constants  and  i  I  ton.  When 
light  rays  shine  through  ellipse  /',  the  equation  for 
their  projection  is 

(t/0,)*’  Mftt,)2  (f/r)’.  (A2) 


When  that  surface  intersects  the  plane  of  the  record 
at  f  //,  the  resulting  curve  is  as  follows: 

(WilM/f/a,)2  =  (tlr)2. 

This  last  equation  can  be  rewritten  as  a  hyperbola 
equation  in  standard  form: 


Comparing  this  equation  for  the  projected  curve 
with  equation  7  gives  the  parameters^  and  by  which 
are  required  for  the  velocity  solution. 


--  HtI a, 

(A4) 

--  «fc/<V 

(AS) 

/  tqure  M  A  hyperbola  pro/ei  tor  tor  measuring  reflection  eelocittei. 


c* 


i » r 


$id«  Vitw 


/  iqure  A  2.  /  he  hyperbola  pro/et  tor  qeometr  v. 
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These  two  equations  show  that  the  design  of  the  pro¬ 
jector  is  based  on  an  appropriate  choice  of  the  height 
H  and  the  etching  of  a  suitable  range  of  ellipses  with 
constants  a,  and  0j.  The  remaining  parameter,  t, 
must  be  measured  during  the  curve-matching 
operation. 

A  calibration  test  on  the  projector  confirmed  that 
there  were  no  significant  differences  between  the  de¬ 
sign  and  the  performance  of  the  projector.  The  reso¬ 
lution  of  the  method  was  estimated  by  using  one-half 
the  difference  in  velocity  of  adjacent  curves.  Accord¬ 


ingly,  for  the  oil  exploration  records,  the  resolution 
was  ±0.1  km/s.  For  the  sub-bottom  profiles,  it  was 
±  0.25  km/s. 

The  velocities  from  sub-bottom  profiles  have  an 
additional  source  of  possible  error  arising  from  their 
dependence  on  ship's  speed.  According  to  the  deriva¬ 
tion  in  the  Methods  section,  the  seismic  velocity  esti¬ 
mate  is  directly  proportional  to  the  ship’s  speed,  and 
therefore,  a  possible  error  of  ±  10%  in  the  ship’s 
speed  causes  a  corresponding  error  in  the  velocity 
determination:  ±0.15  km/s. 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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